H_29_11

Marcus Groll

27 November 2016

Slide with R Output

install.packages("raster", repos = "http://cran.us.r-project.org")
## Installing package into 'C:/Users/Marcus/Documents/R/win-library/3.3'
## (as 'lib' is unspecified)
## package 'raster' successfully unpacked and MD5 sums checked
## Warning: cannot remove prior installation of package 'raster'
## 
## The downloaded binary packages are in
##  C:\Users\Marcus\AppData\Local\Temp\Rtmpk9AeZM\downloaded_packages
library(raster)
## Loading required package: sp
install.packages("cluster", repos = "http://cran.us.r-project.org")
## Installing package into 'C:/Users/Marcus/Documents/R/win-library/3.3'
## (as 'lib' is unspecified)
## package 'cluster' successfully unpacked and MD5 sums checked
## Warning: cannot remove prior installation of package 'cluster'
## 
## The downloaded binary packages are in
##  C:\Users\Marcus\AppData\Local\Temp\Rtmpk9AeZM\downloaded_packages
library(cluster)

NDVI <- function (NIR,red){
  (NIR-red)/(NIR+red)
}

band3 <- raster("F:/Uni/EAGLE/MB1/01_Data/02_Raster/Landsat-5_Ammersee/LT51930272003195MTI01/LT51930272003195MTI01_B3.TIF")

allbands  <- brick("F:/Uni/EAGLE/MB1/01_Data/02_Raster/RapidEye_Ammersee/re20110906_ortho_ls_prj_grto0509_sub.tif")

Rapid_NDVI <- NDVI (allbands[[5]],allbands[[3]])

Slide with Plot

plotRGB (allbands, 3,2,1)

plotRGB (allbands, 3,2,1, stretch ="lin")

plot(Rapid_NDVI)